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We study the issue of the selection of viscous fingering patterns in the limit of small surface 
tension. Through detailed simulations of anisotropic fingering, we demonstrate conclusively that no 
selection independent of the small-scale cutoff (macroscopic selection) occurs in this system. Rather, 
the small-scale cutoff completely controls the pattern, even on short time scales, in accord with the 
theory of microscopic solvability. We demonstrate that ordered patterns are dynamically selected 
only for not too small surface tensions. For extremely small surface tensions, the system exhibits 
chaotic behavior and no regular pattern is realized. 
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There has been a continuing debate regarding the issue of the selection of viscous fingering patterns in the limit 
of small surface tension. On the one hand, the role of surface tension in determining a unique stable steady-state 
finger through the microscopic solvability mechanism jl], ^| has led to the belief that the details of the small-scale 
restabilization control the observed patterns, notably including the selection of a single finger filling half the channel 
from a family of such fingers found initially by Saffman and Taylor ||. Other researchers [|, |(| have proposed 
various selection criteria which do not invoke the presence of a microscopic cut-off, and therefore may be entitled 
macroscopic selection. Here, we perform detailed simulations of anisotropic fingering and demonstrate conclusively 
that no macroscopic selection occurs in this system. 

Viscous fingering, wherein a low viscosity fluid invades a more viscous one in a Hele-Shaw channel geometry, (see Fig. 
|l|) , first described by Saffman and Taylor || , is one of the classic pattern forming systems and serves as the basis for 
much of our intuition regarding pattern formation. Saffman and Taylor showed both theoretically and experimentally 
that this system exhibits a pattern-forming instability, an instability that is restabilized at short length scales by the 
, action of the surface tension between the two immiscible fluids. Detailed studies (IJ fj], |8[ || of steady-state propagating 
finger solutions in this system have revealed the following: At zero surface tension there exists a continuous family of 
solutions, with continuously varying velocity and asymmetry. At any finite value of the surface tension, 7, this family 
is reduced to a discrete infinity of symmetric solutions, of which only the most narrow (and fastest) is linearly stable. 
As 7 is reduced toward zero, all these solutions approach a width that is a fraction A = 1/2 of the channel width. The 
physical interpretation of these results has proven very controversial. The claim has been made that it is physically 
unreasonable for extremely small surface tension to "select" the A = 1/2 finger in times of order 1, from a continuum 
of fingers whose curvatures are never large and should therefore all by essentially unaffected by the surface tension. It 
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FIG. 1: Schematic of the viscous fingering setup. A Hele-Shaw cell is filled with a viscous fluid (here water) and a less viscous 
fluid (here air) is injected uniformly from one end. The result is a finger propagating down the channel (in the direction of 
increasing y). 
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FIG. 2: (a)Tip curvature ko vs. t for various values of isotropic surface tension, 7, starting from identical initial conditions of 
a small symmetric perturbation. (b)Same for case of crystalline anisotropic surface tension with e = 0.4. The number of points 
on the interface, N, was taken between 200 and 800, depending on 7. The local relative error of the Gear's method integrator 
was fixed at 10 _ 6 



is been argued that some other, macroscopic, criteria, must then be responsible for the dynamically chosen A = 1/2 
finger. We will demonstrate in the following that these claims are incorrect, and indeed the microscopic regularization 
provided by the surface tension is dynamically relevant and indeed controls the pattern selection. 

Part of the difficulty in unraveling the puzzle of selection is that both proposed selection paradigms lead to the same 
A = 1/2 finger at small surface tension. It is thus useful to consider a situation where the two mechanisms would lead 
to different results. Such is the case of viscous fingering with a crystalline anisotropy of the surface tension. Here the 
boundary condition for the pressure on the interface, Pint, instead of the standard p- m t = — JK, (for surface tension 7 
and k = —y"(x)/(l + y'(x) 2 ) 3 ^ 2 being the curvature of the interface), takes the form 

Pint = — 7o(l — ecos46*)>t (Crystalline anisotropy) (1) 

where 9 is the angle made by the interface normal to the overall flow direction, y. This form is familiar from studies of 
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FIG. 3: Simulation of Ko(t) using cutoff (B), with 8 = 1 and varying 71. The method and initial condition are the same as in 
Fig. I 
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dendritic crystal growth [g, |HJ and has already been used to model experiments with artificially imposed anisotropy 
Pt [l3[ [TJ] . Study of the steady-state finger with crystalline anisotropy has revealed [jl(| |l5) that for any positive 
e, the unique stable finger is narrowed, with a width of order 7q^ 2 for small 70 (a scaling identical to that of dendritic 
solidification of a solid with crystalline anisotropy). Investigation of the dynamics of such an anisotropic system 
should thus reveal the competition between a macroscopic selection mechanism, if it existed, with its favoring of the 
A = 1/2 finger, and the surface tension which now favors a very narrow finger. We do this by solving the initial 
value problem (for an inviscid pushing fluid) to see when and how the steady-state is achieved. In Fig. ||, we present 
the results of this simulation, graphing the tip curvature, kq, as a function of time for varying 7. The simulations 
were performed via the boundary integral technique described in Refs. |0, |l2|. The interface is parametrized by 9(a) 
where < a < 1 is the relative arclength, and reflection symmetry of the finger is assumed. The initial interface was 
chosen slightly perturbed, 

9{a) = 0.2sin(7ra) . (2) 

In Fig. ^i, the results for isotropic surface tension are presented. We see that kq quickly rises, saturating at a 7 
dependent value slightly below the value of it which characterizes the A = 1/2 Saffman- Taylor finger. (The Saffman- 
Taylor finger is given by the curve y(x) = 2( - 1 ~ A ^ lncos(?r) )• The selection occurs on a time scale essentially 
independent of 7, leading apparent credence to the idea that the selection mechanism is independent of 7. In Fig. ||b, 
however, we present the results of the same simulation, now performed with crystalline anisotropy, with the anisotropy 
parameter fixed at e = 0.4. If the selection was driven by macroscopic effects, the pattern should be insensitive to 
the precise form of the surface energy for short and intermediate times. For such times, the data should, for short 
and intermediate times, recapitulate that of Fig. ||a. Only at long times, times that diverge as 70 — > 0, should the 
narrow finger of the steady-state theory emerge. This is not at all what occurs. Instead, the selection of a narrow 
finger with large tip curvature actually is faster than for the isotropic case, with the selection time actually decreasing 
with decreasing 7. This is also consistent with the experimental finding that imposing anisotropy creates dendritic 
structures on a fast time scale Jll| . 

The generalization for dynamics of the microscopic solvability theory of the steady state that emerges from these 



simulations (see also Ref . 16 ) is that as the regularization embodied in the pressure boundary condition is removed, 
the system tracks some particular time-dependent zero-surface tension solution (r|; but, one cannot determine which 
of those solutions is selected by any macroscopic construction. The existence of multiple trajectories emanating from 
arbitrarily close initial conditions is a consequence of the Hadamard ill-posed nature of the7 = problem Fig. ^> 
shows that the selected dynamics with crystalline anisotropy is well-approximated by the curve Ko(t) ~ C(t* —t)~ 2 until 
such time as it saturates to kq(oo) ~ 7o~ 5 - This behavior is exactly what is expected based on the above principle. As 
pointed out by Howison Jl^j and also by Shraiman and Bensimon [ [l9|| , the zero-surface tension dynamics generically 
exhibits finite-time singularities leading to cusped interfaces with infinite curvature. It was shown that these cusps 
can be characterized as places where the conformal map z = x + iy = F(£) which maps the exterior of the unit circle 
to the actual flow domain becomes singular via having a zero of ^ reach the unit circle from its interior. Assume 
this happens at time t* and, without loss of generality, at zero phase; we will also assume for simplicity that there is 
no tangential component of the singularity velocity, since we are dealing with symmetric initial conditions. (It is easy 
to check that this restriction changes nothing.) Near this time, we have [l!| 

F{e ie ) ~ K (e ie - -) . (3) 



This leads to the parametric curves 



x(0) ~ -K(2-t/t*)6 2 + constant 

y(0) ~ K (0(1 - t/t*) - 9 3 ) (4) 

Hence, the curve approaches the cusp structure \x\ ~ \y\i via the blow-up kq ~ x"(y) ~ C/(t* — t) 2 . Thus, this form 
of anisotropic surface tension has selected the zero-surface tension solution with a finite-time singularity, rather than 
the A = 1/2 finger. Moreover it has done so on a time scale independent of 7. Clearly which pattern is dynamically 
selected (even on short time scales) is controlled by the surface tension, despite its smallness. 

We have seen that we can, by picking a suitable surface tension, cause either the A = 1/2 or the near-cusp solution 
(with tip curvature proportional to 7 -1 ^ 2 , although 7K at the tip is still small) to be selected. In fact, by a suitable 
form of surface tension anistropy, for example, 

y"(x) 

Pint =71 T, ; 7, xon3 (Mesoscopic anisotropy) (5) 

[1 + + y'(xy)2 
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FIG. 4: Simulation of Ko(t) for isotropic surface tension 7 = 0.0002 and the same initial condition as in Fig. ^. Here N = 1000 
and the relative error 10 -9 . 



we can select any A(< 1/2) we choose. In Fig. ||, we present our results for this mesoscopic anisotropy, showing that 
now the system chooses, in finite time, a finger with tip curvature which approaches a finite value significantly greater 
than 7r. Again, there is no sign of a macroscopically selected finger of width A equal to 1/2. To understand our result, 
we note that solvability theory [jio) predicts that the width is determined by matching the complex-plane singularity 
in the modified curvature that occurs when y'(x) = iyl + 6 with a zero of the function g(y') = y"(y'(x)). For the 
finger, 

\\ / / W \ 2 \ 

(6) 




vanishes at y' = i(l — A)/A, giving the selected values 

A = 



l + VT+5 ' 

k = ^VTTS(l + y/T+S) . (7) 

This result is in good agreement with our simulations and rules out any macroscopic selection. Rather, the detailed 
form of the microscopic re-stabilization indeed determines the selected pattern. 

However, by no means is this the whole story. We saw in Fig. ||a how reducing 7 for the usual isotropic surface 
tension appeared to yield a limit curve which asymptotes to the selected Ko = 7r corresponding to A = 1/2. Let us 
examine what happens for an even smaller 7 with exactly the same initial conditions. We see in Fig. |4| that k starts 
out as before, but then turns down, going negative, due to a splitting of the tip |(| [H). How are we to understand 



this, given the demonstration 122, 23, B4L E5| of the linear stability of the selected finger for all surface tension? The 



explanation is that the size of the basin of attraction of the A = 1/2 steady state finger is 7 dependent, becoming 



smaller (exponentially quickly) as 7 — > 123, 23 



Thus, the complaint that extremely small surface tension is incapable of selecting the finger width is perfectly 
correct - at such small surface tensions, indeed, no steady-state finger is selected. Rather, it is likely the system is 
chaotic, generating a random succession of tip-splittings similar to that seen in the Diffusion-Limited- Aggregation 
model pq | of Laplacian growth, where noise is added explicitly. We can see this by examining the divergence of two 
nearby trajectories, as shown in Fig. |^. Two runs with very slightly different initial conditions were performed, the 
first as in Eq. || and the second with the size of the initial perturbation 0.001% larger. We find that initially the 
difference in k,q (t) slowly increases in a linear fashion, but eventually rises exponentially, increasing by over two orders 
of magnitude. Presumably, as 7 is reduced to zero, the maximal Lyapunov exponent diverges, giving rise to the 
Hadamard ill-posedness of the zero 7 problem. Selection of the Awl/2 finger is thus seen as a kind of intermediate 
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FIG. 5: Difference of Kq{€) for two runs with isotropic surface tension 7 = 0.0002 and slightly different initial conditions. Here, 
N = 200 and the relative error 10" 8 . 

asymptotics, realized for small but not too small 7. Those searching for a macroscopic selection principle to account 
for selection of A = 1/2 at extremely small surface tension are then searching for the solution to a nonexistent dilemna. 
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